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Abstract 

The production of energetic particles in the universe remains one of 
the great mysteries of modern science. The mechanisms of acceleration 
in astrophysical sources and the details about the propagation through 
the galactic and extragalactic media are still to be defined. In recent 
years, the cosmic ray flux has been measured with high precision in 
the energy range from 10^*^ to eV by several experiments using 

different techniques. In some energy ranges, it has been possible to 
determine the flux of individual elements (hydrogen to iron nuclei). 

This paper explores an astrophysical scenario in which only our Galaxy 
and the radio galaxy Cen A produce all particles measured on Earth in 
the energy range from 10^° to eV. Data from AMS-02, CREAM, 
KASCADE, KASCADE-Grande and the Pierre Auger Observatories 
are considered. The model developed here is able to describe the total 
and individual particle flux of all experiments considered. It is shown 
that the theory used here is able to describe the smooth transition 
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from space-based to ground-based measurements. The flux of each 
element as determined by KASCADE and KASCADE-Grande and 
the mass sensitivity parameter A^ax measured by the Pierre Auger 
Observatory above 10^® eV are also explored within the framework of 
the model. The transition from 10^® to 10^® eV is carefully analyzed. 

It is shown that the data measured in this energy range suggest the 
existence of an extra component of cosmic rays yet to be understood. 

Keywords: cosmic rays, galaxies: active - galaxies, galaxies: jets - galaxies: starburst 


1 Introduction 

The currently accepted view is that the cosmic rays are produced in ac¬ 
tive astrophysical objects: supernovae, gamma-ray bursts, active stars in 
binary systems, pulsars, active galactic nuclei, quasars, radio galaxies, and 
large-scale structure shocks. Possible sources of cosmic rays in the Galaxy 
include supernova explosions, pulsars and the Galactic nucleus, which con¬ 
tains a super-massive black hole. Traditional models of stochastic acceler¬ 
ation assume the interaction of particles with magnetic fields, according to 
the Fermi mechanisms [T], where the particles would be accelerated in colli¬ 
sions with magnetic moving clouds. In principle, this same mechanism, with 
few changes, can accelerate particles in shock waves in supernovae, gamma- 
ray bursts, Wolf-Rayet star winds, active galactic nuclei, radio galaxies and 
other sites. This paper explores a model originally proposed in reference |2] 
in which a combination of four main components is used to explain the cos¬ 
mic ray spectrum. Figure shows a schematic picture of the main features of 
the model. The details of the model are discussed in the next section. It has 
been previously shown in reference [3] that this model is able to describe the 
total flux of cosmic rays up to 10^° eV using one extra parameter: an energy 
shift factor of 2800 caused by re-acceleration of galactic seed particles in the 
jets. 

The calculations presented here extend the validity of the model in two 
ways. Firstly, an analysis is presented of the cosmic ray flux with energies 
between 10^° and 10^® eV. The recently published data from AMS-02 |1] and 
GREAM [Hin] are used. These experiments are able to discriminate with high 
precision the individual elements of the cosmic ray composition, measuring 
the flux of each particle type. Given the high statistics achieved by these 
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space and balloon-borne experiments, the data constrain severely the contri- 
bntion of each particle to the total flnx. It is shown here for the hrst time 
how this model describes very well the energy range from 10^° to 10^® eV, 
which inclndes the transition regime of space and balloon experiments to 
gronnd-based observatories. This hrst analysis represents an extension in 
the energy range for which this model is able to describe the measnred total 
hnx of cosmic rays. 

In a second analysis, we show the agreement of the model concerning the 
hnx of each element: hydrogen to iron nnclei. For energies below 10^^ eV the 
analysis can be done taking into acconnt the individnal hnx measnred by 
AMS and CREAM for each element. For energies above 10^® eV the analysis 
is done throngh considering indirect composition measnrements. The data 
from KASCADE and KASCADE-Grande experiments are compared to the 
model predictions hxed by the space and balloon experiments. It is shown 
that a continnation extension of the AMS and CREAM (10^° to 10^^ eV) 
to KASCADE and KASCADE-Grande (10^^ to 10^^ eV) data is very hard 
to achieve dne to the high hnx of heavy elements measnred by KASCADE- 
Grande in the energy range from 10^^ to 10^^ eV. The model prediction is 
also compared to the evolntion of the mean depth of the shower maximnm 
(A^max) measnred by the Pierre Anger Observatory [7]. 

Section [previews the original model and its tests. Section shows the hrst 
analysis in which the model is extended to low energies and section shows 
the analysis regarding the cosmic ray composition. Section [^snmmarizes the 
main conclnsions of the paper. 

2 The original model 

Fignre [T] shows the main featnres of the original model which was proposed 
to explain the observed featnres of cosmic rays in the energy range from 10^^ 
to 10^® eV. The model is based on three (1, 2 and 3) galactic and one (4) 
extragalactic component. The phenomena contribnting to the acceleration 
of particles are: a) snpernova explosions into the interstellar medinm, b) 
snpernova explosions into the stellar wind, and c) powerfnl jets of radio¬ 
galaxies. In snmmary, snpernova explosions generate the galactic cosmic 
rays np to eV and radio-galaxies jets re-accelerate galactic cosmic 

rays to the highest energies 10^®“^° eV. In this scenario, the main sonrce of 
extragalactic cosmic rays is Gen A. 
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Label 1 in figure is the resulting energy spectrum outcome of supernova 
explosions into the interstellar medium. The maximum energy that a cosmic 
ray can be accelerated in a supernova shock taking into account the Sedov 
expansion into the interstellar medium was calculated in reference [ 8 ]. The 
produced spectrum has index value proposed to be around —2.75 and an 
exponential cutoff. This component can be written as: 

{dN/dE)i = Til • ^ 

where Ai is the normalization of the flux and jg cutoff energy. 

The cutoff energy is predicted to be proportional to charge oc Z). 

Label 2 in figure is the resulting energy spectrum outcome of supernova 
explosions into the stellar wind, like a Wolf-Rayet star explosion. The pro¬ 
duced spectrum has an index around —2.67 for energies smaller than 
and index around —3.07 for energies smaller than which determines 

the exponential cutoff of the flux. and are both predicted to 

be proportional to charge. The existence of two regimes is due to the de¬ 
pendence of the acceleration efficiency to the particle drift gain [2]. This 
component can be written as: 


{dN/dE)2 


A 2 ■ if E < 

B 2 ■ ■ exp ifE> 


( 2 ) 


where A 2 and B 2 are normalization of the flux. 

Label 3 in figure is an extra component resulting from the outcome 
of a supernova explosions into the stellar wind. In the final stage of the 
very massive stars there is a connection between rotation and magnetic field. 
This magneto-rotational mechanism for massive stars explosions was first 
proposed by Bisnovatyi-Kogan 0 [TO] and seems consistent with the en¬ 
ergy/charge ratio for the heavy elements [TTl [T21 [131 IE!- This connection 
produces a polar cap component relevant in the region where the radial field 
Br ~ 1/r^ dominates. The energy spectrum index was predicted to be around 
—2.33 with a sharp cutoff at ■. 


{dN/dE)s= yl3-E-2-33 ifE<E^'^^°^^ (3) 

where A 3 is the normalization of the flux and is the cutoff energy. 

The original model predicts . 
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This model proposal has used the concept, that transport of cosmic rays 
is governed by Kolmogorov turbulence, and that the secondary particles are 
produced in interactions near the source [TTl [T5] . 

Label 4 in hgure [T] is an extragalactic component proposed to explain 
the highest energy range. Radio Galaxies such as the Fanaroff-Riley class II 
have hot spots at the end of linear radio features, which are considered to be 
highly collimated plasma jets. The evolution of these powerful radio galaxies 
can explain the spectrum to energies above 10^® eV [ISIIIT]. The predicted 
index of the generated energy spectrum is approximately —2. 

It has been shown that the same mechanism rescaled in energy by a factor 
of 2800 can accelerate particles up to 10^° eV |3]. The argument was based on 
the re-acceleration of the original galactic seeds in the jets of radio galaxies. 
Interpretation of observations to derive the central Lorentz factor required in 
the relativistic jets emanating from near super-massive black holes in Active 
Galactic Nuclei (AGN) suggest values of up to = 100 |THl|19]. As Gallant 
& Achterberg [20] as well as |2I] have shown, the acceleration of particles in 
relativistic shocks, clearly possible in AGN jets up to maximally the Lorentz 
factor of the jet itself, gives an increase in energy/momentum by in a 
single hrst step, and for all subsequent steps considerably less. So we use 
here what could be called the “single kick approximation”, namely only that 
single hrst step. Observations suggest that jets are energized intermittently 
(see, e.g., the radio galaxy Her A, [22] )• Such extreme Lorentz factors may 
be possible in the “working surface” of a freshly energized jet. 

3 Comparison of the model to measured en¬ 
ergy spectra of elements 

In previous publications, the predictions of the original model [2] and the 
predictions of its extrapolation to the highest energies [3] were compared 
to the total hux of cosmic ray particles measured by several experiments. 
These comparisons have been able to show the general validity of the model. 
Nevertheless they have not been able to remove the intrinsic degeneracy 
of the model concerning the abundance of each element. If only the sum 
of all elements is verihed, several predictions of the model, for instance, 
rigidity dependencies cannot be tested. Besides that, the large number of free 
parameters to £t the total flux reduces the significance of the final results. 
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In this paper, the tests are done using the most up-to-date spectrum of each 
element as measured by several experiments. The intention of the analysis 
presented in this section is to remove the freedom in describing the total flux 
as presented in the previous studies. 

The original model flux is dominated by component 2 see hgure Com¬ 
ponent 3 causes a break in spectrum which could be used to discriminate the 
model explored here from other traditional models, i.e., Peters cycles [2^12^. 
The step in the energy spectrum caused by components 1 and 2 requires a 
very precise measurement of the energy spectrum of each element in order 
to be tested. 

Data from AMS-02 |1] and CREAM [51E] have been used to validate the 
model in the energy range from 10^° < E < 10^^ eV. Both experiments pub¬ 
lished the H and He flux as shown in ????. The data from both experiments 
offers a very hard constraint to the normalization constants (Ai) for H {A^) 
and He elements. Once these parameters are set to the AMS-02 and 

CREAM data the relative contribution of these elements to the total flux are 
kept constant in the entire energy range studied in this paper (10^° < E < 
10^°eV). Propagation effects, i.e., photo-nuclear disintegration might change 
the relative contribution of each element arriving on Earth at the highest 
energies E > 10^”^ eV. 

The energy cutoff step at cannot be determined due to the lack 

of data in the energy range between and 10^® eV. However, the data mea¬ 

sured by KASCADE [25l [26] with energy above 10^^ eV can be seen together 
with the AMS-02 and CREAM data for H in hgure The hgure the agree¬ 
ment in the relative hux of H as measured by the three experiments is clear 
from the hgure. The agreement of the model to the data is also remarkable. 
The data from KASCADE was used to calculate = 1.96 x 10^^ 

eV which is the energy of the knee of Hydrogen. Using the rigidity depen¬ 
dence of the model = Z x the energy breaks of other 

elements are determined. Figure also shows the data of the KASCADE- 
Grande experiment [271 EH|. The model is able to describe the connection 
between the KASCADE and KASCADE-Grande data which shows a con¬ 
tinuous reduction of the hux up to 3 x 10^® eV. Beyond this energy, the 
KASGADE-Grande data suggest an battening of the proton hux [29]. This 
energy sets the change of predominance from component 2 to component 4 
as show in hgure Again the rigidity dependency model is used to set the 
energy beyond which the extragalactic hux is predominant. 

The rigidity dependency and the relative hux of GNO, NeS and GlMn 
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can be verified and adjusted using the KASCADE and KASCADE-Grande 
data [27]. Figure]^ shows the energy spectrum of intermediate mass particles 
as measured by KASCADE and KASCADE-Grande experiments. The lines 
shown for CNO, NeS and ClMn are the result of the fit of the original model 
to the data. Since the energy breaks have been set by fitting the H spectrum 
and using the rigidity model dependency, the only free parameters in the fit 
are the normalization of each element flux. The agreement of the model to the 
data is very good. The spikes in the flux of each element caused by 
are visible in the sum of all intermediate elements with small amplitude. 
Unfortunately the resolution of the measurement is not enough to test the 
small spikes. 

Finally the model was compared to the iron flux measured by CREAM, 
KASCADE and KASCADE-Grande. Figure]^ shows the fit of the model to 
the data considering two approaches. First the rigidity dependency of the 
energy cutoff was kept = 26 x as shown by 

the dashed line. The label “This model - Rigidity Dependency” is used to 
identify the hypothesis in which = Z x for all elements. 

It is clear that in this case the model does not describe the KASCADE 
and KASCADE-Grande data. If the rigidity dependency is not kept and 
j^cutoff-Fe allowed to vary, the model describes fairly well the data as shown 
by the full line in figure [s In this case the fit leads to = 2.54 x 10^^ 

eV. The label “This model - Fe excess” is used to identify the hypothesis in 
which = Z X for all elements except Fe. 

This analysis illustrates two possibilities to this model and at some extent 
to any rigidity dependent model conceived to describe the energy spectrum 
of cosmic rays with energy between 10^® and 10^® eV. The energy spectrum 
of iron nuclei measured by KASCADE and KASCADE-Grande seems to 
require an extra flux of heavy particles for energies between 10^^ and 10^® eV. 
This extra flux can be provided if the rigidity dependency of the knees is not 
maintained or if an extra flux of iron from a yet unknown source is produced. 

The prediction of the model as fitted to the energy spectra of the ele¬ 
ments shown in figure ?????? was summed in order to obtain the total flux 
of particles. Figure]^ shows the energy spectrum of all particles as measured 
by KASCADE |25|, KASCADE-Grande [29| and The Pierre Auger Obser¬ 
vatory [30]. The spectra of each element as measured by CREAM (H, He 
and Fe) and AMS-02 (H and He) are also shown. Two possibilities for the 
total flux predicted by the model are shown. The full black line takes into 
account the iron flux that fits the KASCADE-Grande data better, which 
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does not obey the rigidity dependency of the knee. The dashed black line 
takes into account the iron flux that hts the KASCADE-Grande data worse 
which retained the rigidity dependency of the knee. 

4 Comparison of the model to the depth of 
shower maximum 

For energies above 10^® eV, the most reliable composition parameter is the 
depth of the shower maximum (Xmax)- The results published by the Pierre 
Auger Collaboration [21] concerning the evolution of the X^ax with energy 
is independent of shower simulation and detector efficiencies, therefore this 
datum is used for comparison to the model predictions. The KASCADE and 
KASCADE-Grande data can also be transformed to an equivalent (Xmax)- 
Using the parametrization of (X^ax) as a function of energy and mass pub¬ 
lished in reference [32| and the flux of element groups measured by the exper¬ 
iments [25] EE] [33| it is possible to calculate the equivalent (Xmax) of the mea¬ 
sured abundance. Figure shows the comparison of the model as optimized 
in the previous section to the transformed KASCADE data, transformed 
KASCADE-Grande data, and Auger data. The energy spectra predicted by 
the model was transformed into (X max ) measurements using the same pro¬ 
cedure adopted for the KASCADE and KASCADE-Grande data, using the 
parametrization of (X ma x) as a function of energy and mass published in 
reference [32] for the Sibyll hadronic interaction model [3^ . 

Both versions of the model with and without an extra iron component 
describe the data qualitatively well. The Fe excess artihcially introduced 
around 10^^ eV produces a decrease of the (Xmax) and makes the transi¬ 
tion between the transformed KASGADE-Grande data and Auger data more 
abrupt. 

5 Final Remarks 

The data analyzed here conhrm the validity of the model in which only the 
Galaxy and GenA could produce all cosmic rays measure on Earth with 
energy above 10^° eV. For the hrst time, the data from AMS-02 and GREAM 
were used to £x the relative contribution of individual elements. 


The analysis presented here tries to describe the iron nnclei flnx recon- 
strncted by the KASCADE-Grande experiment within the framework of the 
original model [2]. The calcnlations shown in the previons section snggest 
an extra flnx of iron nnclei which does not obey the rigidity dependence 
hypothesis. The KASCADE-Grande Collaboration reported a snppression 
in the flnx of heavy elements at 8 x 10^® eV |33]- However, according to 
this pnblication the snppression is less signihcant in the all particle spec¬ 
trum [33]. At the same time, it is clear that the steepening of the spectrum 
after the knee is more severe for the light component rather than for the 
heavy component. The calculation done here suggests that this energy range 
1016.5 < eV might contain an extra flux of a heavy element. This 

idea has been advocated by Hillas [35|. However, in his proposal. Galactic 
magnetars are considered as possible sources for the extra iron flux. An¬ 
other similar study was presented in reference |36] in which the need of an 
additional extragalactic component is pointed out. 

Depending on how many supernova explosions contribute cosmic ray par¬ 
ticles near the knee region and beyond, the possibility cannot be excluded 
with certainty, that just relatively few explosions contribute, and then one is 
the strongest. If this strongest one derives from a more powerful supernova, 
such as in the hyper-nova model discussed in reference [371 EH] then the cos¬ 
mic ray spectrum near and beyond the knee region could be modihed. This 
could entail, that the corresponding knee energy, the spectral index below 
the knee, the polar cap component, and the component beyond the knee 
could all be changed, especially by being shifted to higher energies, and that 
the deeper layers of the stars might be exposed, allowing Fe to become much 
stronger for this one star. 

The reconstructed particle fluxes published by KASGADE and KASGADE- 
Grande were converted to an equivalent (Xmax) as shown together with the 
measured (X ^a x) published by the Pierre Auger Gollaboration. Despite the 
large error bars of the converted KASGADE-Grande data for energies above 
10^^'^,the apparent continuity of the KASGADE-Grande to Auger data is 
remarkable. The hypothesis with and without an extra iron nuclei flux for 
the model proposed here are compared to the data. The data cannot to 
discriminate between the two hypothesis. Here we see a possible trade-off, a 
special Fe-component, well supported by data, and a normal Fe-component 
(i.e. as in the rigidity dependent model), and a much enhanced high energy 
H-component. New composition data from the LOFAR, TA and Auger ex¬ 
periments is expected to the published in the near future for the energy range 


9 


10^^ to 10^® eV. These new data might be able to discriminate between the 
two hypotheses. 
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Figure 1: Schematic representation of the energy spectrum as predicted by 
the original model [2]. Four curves are due to the four components: 1) 
supernova explosions in the interstellar medium (Sedov phase), 2) supernova 
explosion into the stellar winds (Wolf-Rayet), 3) polar cap component of 
supernova explosion and 4) extragalactic contribution. is the cutoff 

energy for component 1, marks a step in the flux of component 2 due 

to a regime transition in the drift gain, and is the cutoff energy for 

component 2. Figure adapted from reference j2]. 
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Figure 2: Hydrogen nuclei flux as a function of energy. The CREAM EE], 
AMS-02 |1], KASCADE |25l [26] and KASCADE-Grande [23, |2H] data are 
shown together with the prediction of the model considered here. The model 
was fit to the data and the extrapolation to the highest energies was done 
following reference |3| 


16 




















Figure 3: Helium nuclei flux as a function of energy. The CREAM EE] 
and AMS-02 |3] data are shown together with the prediction of the model 
considered here. The model was £t to the data and the extrapolation to 
the highest energies was done following reference |3]. All energy breaks and 
cutoffs follow the rigidity dependency after the hydrogen £t shown in £gure|^ 
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Figure 4: Intermediate nuclei (He, CNO, NeS and ClMn) flux as a function of 
energy. The CREAM [5l|6], KASCADE [25l|26] and KASCADE-Grande [271 
I2H] data is shown together with the prediction of the model considered here. 
The model was fit to the data and the extrapolation to the highest energies 
was done following reference [3]. All energy breaks and cutoffs follow the 
rigidity dependency after the hydrogen fit shown in figure 
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Figure 5: Iron nuclei flux as a function of energy. The CREAM O E], 
KASCADE [25l[26] and KASCADE-Grande [271 EH] data are shown together 
with the two hypotheses based on the model considered here. “This model - 
Rigidity Dependency” is used to identify the hypothesis in which = 

Z X for all elements. “This model - Fe excess” is used to identify 

the hypothesis in which = Z x for all elements except Fe. 

The model was fit to the data and the extrapolation to the highest energies 
was done following reference [3]. 
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Figure 6: All particle flux as a function of energy. The CREAM [3 [6], AMS- 
02 [1], KASCADE |25l |26] and KASCADE-Grande |23 |28] and Auger [30] 
data is shown together with the prediction of the model considered here. 
The model was fit to the data and the extrapolation to the highest energies 
was done following reference [3] . “This model - Rigidity Dependency” is used 
to identify the hypothesis in which = Z x for all ele¬ 

ments. “This model - Fe excess” is used to identify the hypothesis in which 
j^cutoff-e — 2 ^ j^cutoff-H elements except Fe. The model was fit to 

the data and the extrapolation to the highest energies was done following 
reference [3|- 
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Figure 7: Mean depth of shower maximum ((Xmax)) as a function of energy. 
KASCADE [2S1 EH] and KASCADE-Grande |2Z1 12H] data were converted 
to (Xmax) using reference [22] • The Auger data are also shown [2T]. “This 
model - Rigidity Dependency” is used to identify the hypothesis in which 
j^cutoff-e _ X for all elements. “This model - Fe excess” is used 

to identify the hypothesis in which = Z x for all elements 

except Fe. 
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